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ABSTRACT 

The general theory of MST radar observations of gravity wave spectra is developed in this 
paper. This effort extends the previous results to include anisotropy and Doppler effects for the 
spectra, as well as the consequences for the multibeam configuration. The relationships between 
the observed one- or two-dimensional spectra for the line-of-sight velocity in the gravity wave 
spectra are derived. Expressions for cross spectra, as well as covariances between velocities 
observed on different beams, are computed. Using these results, studies are carried out to show 
how model parameters for gravity wave spectra can be estimated from the observed quantities. 
Model parameters include the variance, power-law indices, anisotropy parameters, Doppler 
parameters, mean scale sizes, etc. Cases with different numbers of beams will be investigated. 

INTRODUCTION 

Recently, MST radars have been used by many investigators to study the spectra of 
velocity fluctuations in the middle atmosphere [Balsley and Carter, 1982; Meek et al., 1985; 
Scheffler and Liu, 1985; Scheffler and Liu, 1986; Smith et al., 1985; Franke et al, 1988, 
VanZandt, 1985]. Attempts have been made by these authors to interpret the observed spectra in 
terms of internal gravity waves. Based on the Garrett-Munk model spectrum for ocean waves, 
VanZandt [1985] proposed a similar model spectrum for atmospheric waves. Features of the 
model spectrum have been compared to observed spectra [Smith et al., 1985; Franke et al., 1988]. 

In this paper, we generalize the model spectrum for the atmosphere to include different 
types of asymmetry in the model. We then derive relations between radar observed quantities and 
the model parameters. Based on these relations, we propose several observational configurations 
and data processing schemes in order to determine, experimentally, some of the parameters in the 
model spectrum. 

MODEL GRAVITY WAVE SPECTRUM 

For a general GM type spectrum for gravity waves the following form is assumed: 

E(<(>,k z .co) = Eo A 2 (k 2 ,o>) B(co) F(<|>,©) (1) 

E 0 is the energy density constant; k z * is the vertical wave number bandwidth; p and t are the 
frequency and wave number power law indices, respectively. The corresponding three algebraic 
model functions are defined as 

A z (k z> 0» = ^l + ^ fe Sv )^Xl)‘ (2a) 


B(o» = 


P~1 

©il-P- ©b^P 


Icol -P for ©i <l©l < ©b, and B(©) = 0 elsewhere (2b) 
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F(<M>) = ^ + -^ (2c) 

1 — — (l-a^) 1 ^cos(<(M|> a ) 1 — — ( 1 — a^) i^cos(<()-<J) a ) 

kol kol 

where a = (1-e 2 ) 1 / 2 . 

The A z function is based upon the familiar GM 75 vertical wave number function 
introduced in the symmetric model with bandwidth k z *. However, the vertical symmetry faction 
S v , has been added to model the upward/downward propagation asymmetry of gravity waves for 
two different cases. The spectral shape is the same for the ± k z (upward/downward) branches, but 
the levels of these two branches are altered so that they are proportional to 1 ± S v . A physically 
realistic spectrum requires that IS V I ^ 1. The extreme example of S v = ± 1 therefore corresponds to 
solely upward/downward propagation. Any case of S v * 0 adds one degree of freedom to the 
model to represent vertical asymmetry. 

The B function retains the same symmetric power law form as previously assumed. In 
comparison to Garrett and Munk, however, a double-sided form is emphasized throughout this 
paper to be consistent with even symmetry of the four-dimensional wave energy spectrum 
necessary to represent a real velocity field in the space-time domain. This physical requirement is 
also the reason for introducing the sign of O) dependence to the A and F functions in order to model 
the vertical and azimuthal asymmetries. Complete separability is therefore not possible in the 
asymmetric spectrum, but 4>-k z separation is still assumed. If so desired, however, the model 
presented here is simple enough to be represented as the sum of four completely separable products 
when the ± branches of k z , and co are treated individually. 

The F function has been introduced to model the azimuthal dependence of the gravity wave 
spectrum. This function is essentially the sum of two ellipses in the azimuthal plane whose pairs 
of foci lie on opposite sides of the origin but along the same axis inclined by the angle <J> a from the 
x-axis. The eccentricity of both ellipses is e, which is determined by the azimuthal shape factor a. 
The azimuthal symmetry factor Sh controls the relative sizes of the two ellipses in varying degrees, 
from waves propagating with bilateral symmetry, Sh = 0, to those propagating only toward one 
axial direction, Sh = 1- The minimum value for the shape faction is a = 0, which represents 
propagation only in the vertical plane inclined by fa from the x-axis. Nonzero shape factors lead to 
general three-dimensional spectra. The case for a maximum value of a = 1 corresponds to 
azimuthal isotropy for which F(<|>,a>) = 1. 

The following normalizations apply to the general model for the functions A, B, and F: 

oo ©o (Oh 0 

J J A z (k z ,co) dk z = J J B«o) dto= J B(co) dco=^- J F(<M) d* - 1 (3) 

— CDj 2n 

Finally, the spectrum becomes the symmetric Boussinesq model spectrum for the special 
case of S v = 0 and a = 1. 

GRAVITY WAVE SPECTRA AND VARIANCES OBSERVED BY RADAR 

We shall proceed to derive the relation which involve the observed frequency spectrum and 
velocitv variance from a given radar beam and the wave spectrum E(Tc,co). For a positioned radar 
beam, r = r B, with range r and beam direction 

A AAA 

B(0b»<1>b) = cos 4 >b sin 0 b X + sin <t>B sin 0 b y + cos 0 b z (4) 


By closely following the procedure in Scheffler and Liu [1985], the observed velocity 
fluctuation along the beam direction is found to be 
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Vob(r.t) = V(rB,t) • B (5) 

where V is the velocity fluctuation associated with the waves, including the three components U, 
VandW. 

From (4) and(5), we can derive the frequency spectrum for the observed velocity 
fluctuation: 

Eob (©) = Ev (<o) cos^b sin^B + Ev (to) sin^B sin 2 0 B + E\y (to) cos 2 0b 

+ ReEuv(to)sin2(t»B sin20B+ ReEuw(to) cos<|)b sin 20 b (6) 

Y+ ReEyw (<o) sin<J>B sin 20 b 

where the rectilinear spectral components Eu(to), Ey(to). Ew(<o), Euv(to), and Eyw(to) are the 
autospectra and cross spectra of the velocity components, respectively. They can be related to the 
general model spectrum for the wave, eq (1), through manipulations following the procedure 
outlined in Scheffler and Liu /3/. The effects of Doppler shift due to background wind can also be 
included [Scheffler and Liu, 1986]. 

The general expression for the velocity covariance Vob 2 at a given radar beam can be 
obtained by integration of the frequency autospectrum 

oo 

<Vob 2 > = f Eob(co)dG) (7) 


which yields 


<Vob 2 > « <lP> cos^b sin^B + <v 2> sin^B sin^Os + <W^> cos^Ob 

+ <UV> sin2<(>B sin^B + <UW> cos<t>B sin20B + <VW> sin<|>B sin20B (8) 


where the rectilinear variance components arc 



<UW> = - Eo Sv’ Sh' “ c°s<t>a 

<VW> = Eo Sv’ S h ’ 


( 9 ) 
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where a* 

© = J©B(©)da) (10) 

©j 

PARAMETER ESTIMATION FOR GENERAL MODEL GRAVITY WAVE SPECTRUM 

Wc will now make use of the variance formulae obtained in the preceding section to outline 
a parameter estimation procedure for the general model spectrum. Let us assume that we have 
simultaneous and continuous Doppler radar velocity time series at the same altitude for at least three 
beams pointing in nonparallel directions. Then the frequency spectra and cross spectra can be 
computed from these beams. Analysis of the spectral slopes can determine the frequency power 
index, p. The Brunt Vaisala frequency, wb, might also be estimated from these spectra or by other 
means such as from background wind and temperature profiles. The inertial frequency can be 
computed from the earth's rotational period and latitude of the radar site. 

Using a 5-beam configuration where beam 1 points upward in the vertical direction, beams 
2 and 4 obliquely to the east and west at an angle qg, respectively, and beams 3 and 5 obliquely to 
the north and south at an angle qg, respectively. If we denote the velocity covariances measured 
by the ith beam as <V| 2 >, then the covariances for the wave associated velocity components can be 
computed from the <Vi 2 >: 


<W 2 > = <Vi 2 > 


<U2> = (<V 2 2 > + <V 4 2 > - 2<W2> cos 2 e B )/(2sin 2 e B ) 
<V 2 > = (<V 3 2 > + <V 4 2 > - 2<w2> cos 2 0g)/(2sin 2 e B ) 
<UW> = (<V 2 2 > - <V 4 2 >)/(2sin 2 0g) 

<VW> = (<V 3 2 > - <V 5 2 >)/(2sin 2 0g) 

If a sixth beam is added, say at = 45°, then it is possible to obtain 


<UV> = 


<y 6 2> _ 1 )(<V 2 2 >+ < V 3 2 >)/44^^-^ 

sin^g 


(<V 4 2 >+<V 5 2 >)/4 


(ID 


With these six measured variances in eq. (11), the solution of eq. (9) allows for the 
quantitative estimation of four parameters for the model spectrum: Eq, a, + tan2<|> a , and IS v hl = IS V I 
Sh. Under the set of assumptions made here, this leaves two degrees of freedom since there are 
six velocity variances in all. Due to this redundancy, it is possible to check the consistence of the 
wave spectrum hypothesis by computing some of the parameters in more than one way. 

Let us first examine the case of E 0 . Inspection of (9) shows that it is determined by <W 2 >, 
or <U 2 . and <V 2 >, 


E° = <W2> * or Eo = ^-( <u2> + <V2>) (12) 

The second way to obtain E 0 from velocity variances is from the horizontal components. If 
<W 2 > is obtained directly from a vertical radar measurement it may be contaminated with a 
horizontal velocity component depending on the possible sloped layer and the precise nature of the 
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scattering mechanism. On the other hand, the value of E 0 obtained by decomposition of oblique 
velocity variance measurements can be compared to the one found by <W^> to check the effect of 
the contamination. Greater consistency could be achieved for the parameter estimation scheme 
applied here if <W^> were also determined by decomposition from an oblique beam configuration. 
The scattering mechanisms for multiple beams at a fixed oblique zenith angle are much likelier to be 
caused by the same process than the vertical beam. This is especially the case of 0g greater than 
about 2.5° where turbulent scatter dominates. A way to obtain <W 2 > from oblique measurements 
is to add the 7th beam at a different beam angle Og* with (fo, say, at 135°. Then 


A sinks’ - B sin 2 9 B 
sin(0 B ' - 0 B ) sin(0 B ' + 0 B ) 


(13) 


where 

A=j [<V 2 2 > + <V 3 2 > +V 4 2 > + <V 5 2 >] 

B = j [<V 6 2 > + <V 7 2>] [<V 3 2> - <V 5 2>] 


Even if the observed velocity variances are obtained under identical scattering mechanisms, 
a discrepancyjbetween Eq computed from (12) could still arise if the mesoscale atmospheric energy 
spectrum E( k ,to) is not entirely due to waves. The presence of two-dimensional atmospheric 
turbulence could contribute to a significant increase in die horizontal variance components without 
having an effect on the vertical variance which would still be due mainly to wave fluctuations. 
This provides a way to examine the contribution of 2-D turbulence to the overall observed velocity 
fluctuations. 

The azimuth width factor, a, can be solved from (9). The most straightforward solution is: 


1-a f«U>-<V 2 >) 2 + (<UV>) 2 ] 1 / 2 

T+a <U 2 xV*> 


Note that for the case of an azimuthally isotropic spectrum, such that a = 1 and so <U 2 > - <V 2 > = 
<UV> = <UW> = <VW> = 0. Then the remaining parameters, <|> a and S v h cannot be determined 
by velocity variances alone. However, if 0 £ a £ 1, then we can continue this approach to solve 
these two parameters also. 

The solution to <|> a is quite straightforward. 


2<UV> . A <VW> 

tanl<t>a = <u^>-<v^> OT tan<t,a= HJW7 


(15) 


The two possible solutions to (|) a cannot be distinguished by variances alone without other 
knowledge or an assumption about the sign of S v . A spectrum which consists of mainly 
downward phase propagation will, for example, have a negative S v , and then (9) can be used to 
determine the quadrant of <)> a . The magnitude of the vertical-horizontal propagation factor S v h can 
be found for a * 1 through the vertical fluxes of the zonal and meridional momenta. 
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I S vh | = I sj S h (j^j n gfcj [<uw2> + <V 2>] 1/2 /E 0 (16) 


This completes all the possible ways to estimate the parameters of the general model 
spectrum by using the nine variance components. Although these variances do include the 
momentum fluxes, they do not determine all the parameters in our model. Spatial information must 
also be available in order to determine kj* and S v (or S h ). A vertical or nearly vertical observed 
wave number spectrum can be directly applied to estimate k z * such as was done by Scheffler and 
Liu [1985]. However, a two-dimensional data set is required to completely determine the 
propagation factors S v or Sj,. The energy density content of each of the two branches of an 
observed k z - co velocity spectrum could, for instance, be used to estimate S v . This parameter 
could also be determined from the same data set, but in the time domain by using the phase 
propagation correlation technique. 

CONCLUSIONS 

In this paper we derived the general expressions for various autospectra and cross spectra, 
as well as different types of variances for velocity fluctuations observed by multibeam MST radars. 
Assuming that these fluctuations are due to the presence of gravity waves, the radar-observed 
quantities are related to parameters of the general model gravity wave spectrum. These relations 
are used to design observational configurations such that the model parameters can be determined 
from the observed data. 

The procedure can be summarized as follows. With one beam, pointing vertically, one can 
determine the spectral indices p and t; also E^ and possibly S v by 2-D spectral analysis. With five 
beams, one vertical, two symmetrically off vertical in the E-W directions, and two N-S directions, 
one can then determine the azimuthal asymmetry parameters, a, and the principal direction of the 
propagation ellipse, With two additional beams, one with the same beam zenith angle but with 
different azimuth angle, one with different zenith angle, sufficient redundancy can be provided for 
checking the consistency of the assumption that the observed velocity fluctuations are due to 
gravity waves. 

Finally, it should be mentioned that the major effort in this procedure is based on the 
variances of the observed quantities. This in many ways is more suitable for the purpose in 
determining model parameters than procedures that use spectral information alone. 
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